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Abstract. We prove the existence of a family of slow decay positive solu- 
tions of a supercritical elliptic equation with Hardy potential in R'^ and study 
stability and oscillation properties of these solutions. We also establish the ex- 
istence of a continuum of stable slow decay positive solutions for the relevant 
exterior Dirichlet problem. 



1. Introduction. 

Our starting point is a superlinear elliptic problem in the entire space 
(1.1) -Au^uP, u>0 inR^, 

where p > 1 and > 3. By ps in what follows we denote the critical 

Sobolev exponent. It is well-known that for p < ps problem (jl.ll) has no positive 
solutions. For finite energy solutions this is an easy consequence of Pohozaev's 
identity. For positive solutions without decay assumptions at infinity this is a deep 
result of Gidas and Spruck . For p = ps all positive solutions of (II. ip are given 
up to translations by a one-parameter family 

Wx{\x\) -.^ X^Wi{X\x\) (A>0), 

iV-2 

where W^i (a;) := (l + {N{N—2))^^\x\'^) ^ is a rescaled minimizer of the Sobolev 
inequality. 

For p > Ps the structure of the solution set of (|l.ip is more complex. First we 
note that for all p > problem (jl.ip possesses an explicit singular radial positive 
solution 
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U^{x) := Cp\x\- — , 

Observe that ii p > ps then Uoo G Hl^^CR^) and hence Uoo is a weak solution of 
(II. ip in the entire R^, despite the singularity at the origin. However Uoo is an 
infinite energy solution because of its slow decay at infinity for p > ps- 

The set of all radially symmetric solutions of (|l.ip can be analyzed through 
phase plane analysis after applying Fowler's transformation, cf. ]1M p. 50-55]. In 
particular, ii p > ps then (jl.ip admits a radial positive solution J7i(|a;|) such that 
Ui{0) = 1. It is known that C/i(|a;|) is monotone decreasing and 

lim ^ ^ C„ 

\x\ P-i 
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however Ui has no expHcit representation in terms of elementary functions. Taking 
into account the scaling invariance one concludes that rescalings of Ui are also 
solutions of (jl.ip . so that (jl.ip possess a one-parameter continuum of radial positive 
solutions 

(1.2) Uxi\x\)=X^Ui{X\x\) (A>0). 

One can show that the singular solution Uao is the limit of the family {U\), in the 
sense that for any x ^ holds 

lim Ux{\x\) = Uoo{\x\). 

A— ^CXD 

In addition, it is known that given < Ai < A2 < 00, solutions U\-^{r) and U\.^{r) 
in the range ps < P < Pjl intersect each other infinitely many times as r — > 00 , 
while for p > pjL the solutions are strictly ordered, that is U\^{r) < Ux^{r) for all 
r > 0. Here 

^-^^^ ifiV>10, 



PJL := { N -A-2^/N~l" 

00 if TV < 10, 

is the Joseph-Lundgren stability exponent, introduced in [TD]. The exponent pj^ 
controls various oscillation and stability properties of solutions U\ , which are par- 
ticularly important in the study of the time-dependent parabolic version of (jl.ip . 
see [51 [T7] or [T3J p. 50-55] for a discussion. 

We are interested in a perturbation of p. ID by the Hardy inverse square poten- 
tial, that is the equation 

(1.3) - Au+ -^u^uP, u>0 inM^\{0}, 

where /i > —Ch and Ch ■— — 4-^^ is the Hardy critical constant, i.e. the optimal 
constant in the Hardy inequality 



(1.4) / \V^\^dx>CH [ j^dx y^eC^iR'^). 

Hardy potential provides an important example of a long range potential, that is 
a potential which modifies asymptotic decay rate of solutions at infinity and their 
behavior at the origin, see e.g. OH]- 

For p PS a. Pohozaev-type identity shows that similarly to (jl.lD , equation (jl.3p 
has no finite energy solutions [16]. For p = ps equation ([LS]) admits an explicit one- 
parameter family of finite energy radial solutions, cf. [21116). However, the structure 
of positive solutions of (jl.3p in the critical regime p = ps is not fully understood. 
It is known that for large values of /i > equation (|1.3p admits nonradial solutions 
which are distinct modulo rescalings from the radial solutions [31 [TB] . See [5] for 
recent results and discussion of open questions in this direction. 

In the present work we consider equation (|1.3p in the supercritical regime p > ^5. 
In the next section we setup the problem and discuss basic properties of the explicit 
singular solution similar to Uoo- In Section [3l for optimal ranges of p and /z we 
establish the existence of a one-parameter family (C/a) a>o of infinite energy solutions 
of (|1.3p . which coincides with (|1.2p when — 0. We also discuss stability properties 
of these solutions. The presence of the Hardy potential produces a range of new 
critical exponents related to stability which do not have immediate analogues in 
the unperturbed case of equation (|l.ip . Finally in Section |4l we discuss equation 
(|1.3p in exterior domains. We justify optimality of critical exponents introduced 
in previous sections. Further, under some assumptions on p and ^ we prove the 
existence of a continuum of infinite energy solutions of (|1.3p , which in some sense 
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could be considered as a perturbation of the original family of solutions on M 
but goes beyond spherically symmetric or scaling invariant setting. This partially 
extends some of the recent results in [4] . 



2. Equations with Hardy potential. 
We study the equation 

2 2 

(2.1) _ Au + ^V^M = , w>0 inR^\^!r, 

where K — {0}, or {0} G K and i^T is a connected compact set with the smooth 
boundary dK, p>l, A^>3, i'>0 and := so that v1 is the Hardy 

critical constant in p.4p . By a solution of (j2.ip we understand a classical solution 
u e C^(R^ \ K), with no apriori assumption on the decay of u{x) at infinity. We 
say u is a weak solution of (12. ip in if m e Hl^^iMJ^ ) and 



j Vu ■ Vipdx + {v^ — j j— ^ dx — j vPipdx 



Note that for v < v^, solutions of (|2.ip must have a singularity at the origin (see 
Lemma [23] below) however this singularity might be compatible with the concept 
of a weak solution in R^. 

We say a solution u of (EJ]) in R^ \ AT has finite energy if m S ^^(R^ \ K), the 
completion of C^(R^ \ K) with respect to the norm || V93||i2. We say a solution u 
of (|2.ip is stable in R^ \ X if the formal second variation at u of the energy which 
corresponds to (|2.ip is nonnegative definite, that is 

(2.2) j \\Iip\^ dx+{u'^-ul) j ^dx-p J uP-\'^dx>Q ^ if e {M.^ \K) . 

A solution u > of (|2.ip is called semi- stable in R^ \ K if it is stable in M.^ \Br, 
for some i? > 0. A solution u > of (|2.1I) is called unstable if it is not semistable. 
Note that these definitions do not require u to be a finite energy solution. 

2.1. Explicit radial solution. For v > Q and p > := 1 H set 
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p-l 



and introduce the critical exponent 

p := 



1 + 17:^' '-^^v<v^, 

oo if 1/ > J/:,, 

Clearly p* > ps- A direct computation shows that Uoo is a positive solution of 
(|2.ip for all < p < p*, while for p ^ [p*,^*] the coefficient Cp^^ becomes negative. 
Note that Uoo e Hl^JR^) for p> ps, that is Uac is a weak solution of (l2J|) in R^. 
However Uoo is an infinite energy solution because of its slow decay at infinity. 

The importance of the solution Uoo is due to the fact that it will be used as an 
elementary building block for constructing further solutions of (|2.ip . In order to 
do this it is essential to understand stability properties of Uoo- 

Lemma 2.1. Let p £ {p*,p*) and v > 0. The solution Uoo is stable if and only if 

(2.3) pC^-' < v\ 

while if (|2.3p fails then Uoo is unstable. 
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Proof. The formal second variation of the energy which corresponds to (|2.ip at Uoo 
is given by 

|V^P dx + {y^ - vl ~ pCP-') j ^ dx. 

Thus the assertion follows directly from the fact that vl is the optimal constant in 
the Hardy inequality (|1.4p . 

Taking into account the scaling invariance of Hardy's inequality we also conclude 
that if (|2.3p fails then Uoo must be unstable. □ 

The inequality (|2.3p amounts to a third degree algebraic expression for which 
closed form solutions could be obtained using Cardano's formulae, however the 
explicit expressions for solutions are tedious. Below we present a qualitative analysis 
of (|2.3p . Set s := —j^, so that (|2.3p transforms into 

1— i < (—1^* — < s < mm{— i^* + i^, (J|). 

\s\ 

Define 

9{s) {s + iy,f{s-2). 
Then solving (|2.3p for < p < p* is equivalent to classifying the roots of the 
equation 

(2.4) e{s) = -2v^ {~v^ - V <s < min{-z/* + v, 0}), 
and solving the inequality 

(2.5) 6{s)<-2v^ - J/ < s < min{-j/* + J/, 0}). 

Note that (?(0) — ~2vl and that d has two critical points: a local maximum at 

c) = and a local minimum at Smin '■— ~ '^'^^ with 
-^(2 + iy,)3. Denote 
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Clearly for every v > equation ()2.4p has exactly one root in the interval 
(— I/* — i/, — z^*). To analyze the roots of (j2.5l) in the interval (— i^, ,min{— i/* + v, 0}) 
we distinguish the cases Smin < and Smin > 0. Q 

In the case Smin > (that is 3 < < 10): 

(i) if 1/ > i^* then (|2.4I) has no roots in (— z/,,0) and (|2.5I) holds Vs <E (— — 

(a) if < i> < v^ then (12.41) has exactly one root a- G (— i^*,— i^* + z^) and 
(P3)) holds Vs e (-i^, - ly, a^] U [(t_, -i^* + v). 

In the case Smin < (that is > 10): 

(i) ii V > V then (|2.4p has no roots in (—;/,, 0) so that (12.51) holds Vs € 

(m) if zy, < < 1/ then (j2.4l) has exactly 2 roots (t_ and ct+ in (— zy,,0) and 
— zy, < (T- < s,„i„ < (T_|_ < so that (12. 5p holds Vs e (— z^, — cr#] U [cr_, cr+], 

iiii) if < z^ < I'* then (j2.4l) has exactly 1 root (t_ in (— zy,,0) and cr_ G 
(— z/*, Smin) so that (|2.5p holds Vs e (— z^* — z/, ct#] U [cr_,0). 

In what follows we denote 

2 2 2 

P# = 1 , p-:=l , P+ := 1 , 

cr# (T_ cr+ 



^Note that if wc write p, = v'^ - ul as in llOt then Ji = TT^iN - W)'^{N - 1). 
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and note that 

l< < p# < PS < P- < P+ < P*, 
for all values of > 3 and v > Q when all the exponents are well defined. Then 
the above analysis leads to the following equivalent to (|2.3|) characterization of the 
stability properties of the solution Uoo in terms of the original parameters p and v. 

Lemma 2.2. Let p e and v > Q. 

(a) // 1^* < v < V and N > 11 then the solution Uoo is stable for p £ U 

. and unstable for p € (p^ , p_ ) U (p+ , p* ) • 

(b) If Q < v < v,, and N > 3 or v = v^, and N > 11 then the solution Uoo is 
stable for p S (p*,p#] U [p-,p*) and unstable for p € 

(c) If v > and 2i<N<lQorv>D and N > 11 the solution Uoo is stable 
for p E (p* , and unstable for p £ [p^ , oo) . 

Remark 2.3. In the pure Laplacian case v = v^, one calculates the explicit values 



N + - 1 N - 2y/N - 1 



iV - 4 + - 1 AT - 4 - 2ViV - 1 

here p- is defined only for iV > 11. Thus for the Laplacian the exponent p_ 
coincides with the Joseph-Lundgren stability exponent, see [TO] or [TH p. 50]; while 
the exponent p^ is known to appear in the context of local singularities of solution 
of equations (|l.ip . cf. [TO] Lemma 5]. 

Remark 2.4. If > 10 and v — v then p- — p^ — -^-zjq is the only supercritical 
value of p where Uoo is stable. 

2.2. Slow and fast decay solutions. Clearly, a solution u of (|2.1I) is a positive 
superharmonic of the linear Hardy operator, that is u satisfy the linear inequation 
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(2.6) - Am + , ,* u > in R'^ \ K. 

m 

As a consequence, solutions of (|2.ip with ly^ < are always singular at the origin 
while for > i/^ solutions might vanish at the origin. More precisely, the following 
local decay properties for positive superharmonics of Hardy's operator hold, cf. 

Lemma 2.5. If u > satisfy (j2.6p in a neighborhood of the origin then 

(2.7) liminf-^^^^ > 0, liminf , "^"^^ < oo. 
If u > satisfy (j2.6l) in an exterior domain then 

(2.8) lim inf , "^^^ > 0, lim inf , f''^ < oo. 

Bidaut-Veron and Veron Theorem 3.3] proved that the structure of the so- 
lution set of (12. 1[) in exterior domains which decay at infinity no slower then Uoo is 
essentially determined by the solutions of the following equation 

(2.9) -Asiv-iw + CP-iw = wP, a;>0 \^S^-^. 
on the sphere S^~^. 

Lemma 2.6. ^ Theorem 3.3] Let p^ps- If u > satisfy (EH]) inR'^\K and 

111 Xj 

(2.10) limsup —'^Koo, 
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then either 

u (x\ 

(2.11) lim „ = c (yasi deca^;, 
|a:|-i.oo |a;| ^ 

or there exists a positive solution lo{-) of (|2.9p such that 

(2.12) lim ^^^M^=w(-) (slow decay) 

\x\^oo p-1 

m i/ie C'^{S^^^) topology, for any k eN. 

Remark 2.7. Clearly, Cp^i^ is a constant solution of (|2.9p . For 1 < p < it is 
known (see [6] or [2l Corollary 6.1]) that Cp^i, is the only solution of p.9p provided 
that 

(2.13) {p^l)C^-' <N-1, 

while if (|2.13p fails then problem (|2.9p admits nonconstant solutions, see [21 Corol- 
lary 6.1], [TBJ Theorem 0.5] and [Ql, Theorem 1.3]. Similar result holds for some 
values p > ^3^, see [3]. The complete structure of solution set of (|2.9p is not yet 
fully understood, see [HI H] for some recent results in this direction. 

Remark 2.8. li ly > Q and p < ps then (|2.10p always holds, see [H Remark 3.2]. 



We will classify positive solutions of (|2.ip into fast and s/ow decay solutions ac- 
cording to alternatives (j2.1ip and (|2.12p . Note that for p > ps slow decay solutions 
are always infinite energy solutions, because of the slow decay rate (|2.12p at infinity. 



3. Radial slow decay solutions in 



pN 



Radial positive solutions u(|2:|) > of (|2.ip in M^\{0} correspond to the positive 
solutions U{r) = u{r) of the initial value problem 

(3.1) -u" -^^—^U' +'^-^U ^U-P (r>0), 

which can be studied through the phase plane analysis. 

The existence of a family of regular at the origin slow-decay solutions of p.ip 
in the Laplacian case v = v^, \s well-known and goes back at least to jlOj . 

Theorem 3.1. Let ps < p < p* . Then for any A > equation (j3.1l) admits a 
unique positive solution U\ G C^(0, oo) such that 

(3.2) linriMrl^A, lim^^C,.. 

Moreover, 

(3.3) Ux{r) = X^UiiXr) VA > 0. 
Further, for A £ (0, oo] the following properties hold: 

(i) if pClpl} < then solutions U\ are stable and ordered in the sense that 
< Ai < A2 < oo implies U\^{r) < U\^{r) for every r > and in addition, 

(3.4) lim U^^^^)-^^^^^) > 0; 

r— s-oc r 

(ii) if pC^I^ > then solutions \J\ unstable and oscillate, in the sense that 
< Ai < A2 < 00 implies that Ux^{r) — U\-^{r) changes sign in {R, +00) 
for arbitrary R > 0. 
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The proof of the theorem fohows the exposition in jl4t pp. 50-53] with minor ad- 
justments needed to accommodate v ^ v^,. We present the sketch of the arguments 
for the readers convenience. 

Proof of Theorem \S.1\ Using the transformation 

(3.5) w{t) ^r-l^U{r), t = log(r), 

problem p.ip becomes an autonomous second order differential equation 

(3.6) w" + 2(3w' + wP - jw ^ 0, te R, 
where since ps < p < p* , 

> 0, and 7 = CP-^ ^ 1^^ - (u, - > 0. 

p-1 P' \ p-lj 



Set 



£{w) = £{w, w') :— —\w' 



2' ' 2 p+l 
Then 5 is a Lyapunov function for (j3.6l) and 

Set X := w and y :— w' . Then p.6|) can be written as an autonomous first order 
system 

'x'\ ^ ( y 
.y' \-2py 7x - xP . 



*(a;,y), 



which possesses two equilibria 

(0,0) and (7^,0) 
in the half-space {{x,y) : x > 0} . Denote 

Ao:=V$(0,0)=(« J^), A.:.V$(7,0)=.(_(^«^^^ 

The matrix has two real eigenvalues 

2 

a± := -P ± \//32 +7 = -~v^±v, 

p - 1 

so that q;_ < < a+. The corresponding eigenvectors are (1, «+) and (1, a_), that 
is (0, 0) is a saddle point of the vector field $. The matrix A^, has two eigenvalues 



:=-/3±v//32-(p-l)7, 
the corresponding eigenvectors are {l,a\) and (l,Q!l). Clearly Re(aj_) < 0, so 
(7?^ , 0) is always an attractor. Note also that a*^ is real if and only if pCP~^ < . 

Using the Lyapunov function E one can show that the trajectory tangent at the 
origin to the eigenvector (1, a-|_) is a heteroclinic orbit which connects the equilibria 
(0,0) and (7~,0), see [H P-52]. Moreover, since (0, 0) is a hyperbolic saddle- 
point, the uniqueness of such hetereclinic orbit follows by standard arguments. 
The corresponding solution w(t) exists for alH G M and satisfies 

(3.7) lim u;(i)=0, lim w(t) ^ . 
Moreover, we can assurae that w(i) satisfies the normalization condition 

(3.8) lim ^ = 1. 
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Since p. 61) is autonomous, w{t + 9) is also a solution of (13.61) that corresponds to 
the same heteroclinic orbit, for any e M. Given G M, set A := e^. Then 

U\{r) := r~^w(log(Ar)) = A^?7(Ar), 

and C/a satisfies p.2p in view of (j3.8p and (j3.7p . that is J/a is the required solution of 
(|3.1|) . The uniqueness of C/a follows from the uniqueness of w{t) since p.Sp defines 
a one to one correspondence between solutions of (j3.ip and (j3.6p . 

To understand oscillation and stability properties of U\ note that the eigenvalues 
a5- are real iff 

/3'> (p-l)7, 

which is equivalent to the stability condition (j2.3p . Note that then 

q;_ < al < —Vi, < a'^ < a+. 

If the roots a^. are real then arguments similar to [TH p. 53] show that the trajectory 
w{t) is monotone increasing in t for all i G R. Hence the solutions U\{r) are 
monotone increasing in A. In particular, U\{r) < Uooir) for any A > and solutions 
U\ are ordered. Further, in view of (j2.3D the solution Uoo is stable. Since U\{r) < 
Uoo{r), we obtain 

By Hardy's inequality we conclude that 

j \V^\''dx+{v^-A) J ^dx-p J Ur\\x\)^'>0 

for all (fi G C^(]R^), that is U\ is a stable solution of (|2.ip . In addition, similarly 
to [m Remark 9.4], we conclude that 

lim = a; > -/3, 

which after returning to the original variables and combined with (j3.3p implies 
(El- 

If Qfj- are complex then similarly to 14, p. 52] one can see that the trajectory 
{x{t),y(t)) spirals infinitely many times around the attractor (7,0) which suggests 
that the solutions U\ oscillate in the sense of (ii). The detailed prove of oscillation 
and instability of U\ when a^- are complex is a particular case of a more general 
Theorem 14.31 which will be proved in the next Section. □ 



Remark 3.2. In the subcritical case < p < Ps equation (|3.ip has no positive 
slow decay solution which satisfy p.2p . Indeed, ii p — ps then (3 — 0, Re{a'^) = 

and the stationary point (7p^,0) is a center. One can show that the trajectory 
tangent at the origin to the eigenvector (l,a+) is a homoclinic orbit. This ho- 
moclinic corresponds to an explicit one parameter family of finite energy solutions 
of dlJ]), see [III pp.253-254]. If < p < ps then ^ > 0, Re(ai) > and the 

stationary point (7P^,0) is repelling. Hence a heteroclinic between (7P^,0) and 
(0, 0) originates at (7^^, 0) and converges to (0, 0) tangentially to the eigenvector 
(l,a_). This heteroclinic corresponds to a positive solution of (13. ip which decays 
at infinity as 0{\x\~'^'~'^) and has a singularity at the origin of order 0(|x|~p^). 
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4. Slow decay solutions in exterior domains. 

First we justify that the value of the nonexistence exponent p* is sharp. The 
result, which is first appeared in [31 Remark 3.2], is an immediate consequence of 
Lemma 12.61 

Theorem 4.1. Let p > p* . Then (j2.ip has no slow decay solutions in \ Bji, 
for arbitrary R > 0. 



Proof. Simply note that for p > p* one has Cp^i, < and hence the equation p.9|) 
on the sphere does not have any positive solution. Then the conclusion follows from 
Lemma 12.61 □ 



Remark 4.2. If p > p* then the slow decay rate is incompatible with the upper 
bound (12. 8p of Lemma 12.51 This argument however does not apply when p = p^. 

Next we justify sharpness of the stability and nonoscillation condition (j2.3l) . The 
result below extends oscillation statement of Theorem IS . 1 1 bevond radial setting. See 
also [171 Proposition 3.5] for related results in the pure Laplacian case v = v^. 

Theorem 4.3. Let p > ps, i' > and pC^^^ > . Let U^, > be a subsolution of 
(1^ such that 

(4.1) liminfiM^ >Cp,,. 

M-^ca p-i 

Then [/, is unstable. Further, if u > is a supersolution of (|l.ip then either 
u = ?7*, or (u — [/»)_ ^ in \ I3]i, for arbitrary large R > 0. 

Proof. From (|4.ip we obtain 

pUr\^)>{i^'+e)\x\-^ {\x\>Rel 

for some e > and R^ > R. Assume that [/, is semistable, that is there exists 
R> such that ((^ holds in \ But then we arrive at 

J |V^p dx + {ly^ - J ^dx>{,y^+e) J ^dx>0 Vy. G C^(R^ \ BrJ, 

a contradiction to Hardy's inequality. We conclude that U^, is unstable. 

Further, set h — u — U^, and assume that ft, > in M.^ \ 13 a, for some R > 1. 
Then by convexity and (j4.ip we obtain 

2 2 

-Ah + " . h > uP -UP ^ iU, + h)P -UP 



X 



VCP'^ ;y2 -I- p2 

> PUr'h='-j^h>^=^h (N>i?e). 

It is well-known that such inequation has no positive solutions, cf. [121 Corollary 
3.2]. We conclude that either ft, = 0, or ft, changes sign in \ Br^. □ 

Remark 4.4. The above result does not exclude possibility that u < in an 
exterior domain. The latter is however not possible in the case when both and 
u are slow decay solutions. In particular, since all the solutions U\ satisfy (14. ip . 
the above result includes the oscillation statement (ii) of Theorem 13.11 



Next we show that if the stability assumption (|2.3p holds then slow decay solu- 
tions of (j2.1l) in exterior domains are well ordered in a certain sense. We consider 
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the exterior boundary value problem for (j2.ip 
(4.2) 



uP, u>0 in \ K, 
u — il) on dK, 



here K ^ {0} is a connected compact set with the smooth boundary dK^ and 
G C{dK) is a nonnegative continuous function. 

Theorem 4.5. Let p > ps, ly > and pG^~^ < . Let U^, > be a slow decay 
solution of (j2.ip in M.^ \ K such that for some R > holds 

U4x)<Uo^{x) {\x\>R). 

Given ip G C{dK) such that 

< V(a;) < C/^x) ondK, 

problem (j4.2p admits a slow decay solution ut such that 

Q<Ut < U., in \ K. 

Moreover, 

(4.3) hm ^ , 

Proof. We are going to construct a sub-solution U_ and a super-solution U such 
that 

<U<U <U^ and C/ = 17 = V on 

Then the existence of a solution uf between U_ and U follows via the classical sub 
and super-solution argument, cf. Theorem 38.1]. 



Subsolution U_. Let > he the minimal positive solution to the problem 
(4.4) ~Ah+ * h:=:pUP-^h mR^\K, h = U^-'4> ondK. 



The existence of such a solution is ensured by the Lax-Milgram theorem. Indeed, 
by the assumptions 

(4.5) pur\x) < pUL-'{x) < vCl:,M'^ < 

Hence the corresponding to (|4.4I) quadratic form is coercive on the Sobolev space 
Dq(R^ \ K). Moreover, from Lemma [2.61 we conclude that given a large i? > 
there exists m £ (0, v"^] such that 

pUr\x) > m\x\-^ {\x\>R). 

Then a standard application of the comparison principle for Hardy operators (cf. 
Lemma [2751 and [12j Lemma A. 8]) implies the two-sided bound 

c|x|"'- <h^< Clxl"- {\x\ > R), 
where a*_ is the smallest root of 

-{a + 1^,. - iy){a + 1^^ + I') = pC^^J^ 

and a'_ is the smallest root of 

— {a + v^: — v){a + v^: + i') = m. 

Note that < m < pC^^^ < i^^, so both equations have real roots and 

, / * 2 
— t^* — ly < a_ < a_ < —v^, < -. 

P- 1 
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Set 

U :— U^, — hfj;. 

Then 

lim . 1, 

and by convexity a direct computation shows 

-A^ = UP - pUr^h^ < ([/, - h^)P = IF in \ ii", 
that is U_ is the required sub-solution. In addition, 

lim El^^lfM^= lin, M^^O, 

|x|-s-oo |x|-!-cx) \x\ 

which implies (j4.3p . 

Supersolution U . Let yy^, > be the minimal solution to the problem 

2 2 

-Ar;+ " h ^ UP-h] inR^XK, r^^U^-^p on OK. 

Note that UP~^ < pUS^^. Hence, solution 77^ exist simply because (|4.5p applies. 
Moreover, a comparison argument similar to the ones above shows that 

< 77^ < /i^ in \ K. 

Define 

U :=?/,- 77^. 

Then 

lim = 1, 

and 

- AI7 - UrHU* -Vi>)> {U* - TyvO""' (t^* - '7^) = C^'' in \ K, 

that is J7 is the required super-solution. □ 

The next result shows that under suitable assumptions on the boundary data the 
exterior problem (j4.2p admits a continuum of distinct slow decay positive solution, 
which in a certain sense could be interpreted as a perturbation of the family of slow 
decay solutions (Ux) constructed in Theorem 13. II 

Corollary 4.6. Let p > ps, ly > and pCP^^ < . Then for every ip G C(dK) 
such that 

(4.6) < V(a;) < Uooix) on OK, 

problem (|4.2D admits a continuum of distinct positive slow decay solutions. 

Proof. Consider the family of slow decay solution (J7a)a>o, constructed in Theorem 
13.11 In view of (14.61) there exists A.0 > such that for all A > 

0<Tp{x) < Uxix) < Uao{x) on dK. 

Let Ai e (A^, 00]. In Theorem 14.51 choose :— and note that in view of (|4.3p 
and p.4p the solution U"^^ given by Theorem 14. 5l is distinct with U\ for any A > A^, 
or with C/^^ for any other A2 > A^. In such a way we have obtained a family of 
distinct slow decay solutions (t^^ )ag(a^,oo] ■ D 



12 VITALY MOROZ AND JEAN VAN SCHAFTINGEN 

Remark 4.7. In particular, ii p > ps and pC^l} < then the problem 

2 2 

(4.7) - Au + , u ^vF ixiM.^\K, u^Q on dK, 

admits a continuum of distinct positive slow decay solutions (C^^)Ae(o,oo] ■ This 
partially extends the result in [H Theorem 1]), obtained in the pure Laplacian case 
V = Vi,. Note however that in [3] the existence of a continuum of slow decay solutions 
was proved for the whole range of exponents p > ps, including the most challenging 
unstable regime pCP~l > p^. The techniques in ^ (see also a survey [5]) are 
based on linearization and perturbation arguments combined with a sophisticated 
machinery of harmonic expansions. Such considerations go beyond the scope of the 
present work. 

Remark 4.8. In the pure Laplacian case = it is known that if K is starshaped 
with respect to infinity then (14.71) has no positive solutions in the subcritical range 
I < p < PS, see [ini Theorem 2]. This suggests that the nonuniqueness statement 
of CoroUarv 14.61 can not be extended beyond the supercritical range of exponents. 
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